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FOREWORD 


This  is  the  tenth  report  in  the  Data  Reduction  Program  Documentation 
series.  It  is  dated  according  to  the  date  of  completion  of  the  documentation.  No 
implication  is  made  that  this  program  will  not  subsequently  be  modified,  amended,  or 
superseded;  on  the  contrary,  the  history  of  radar  data  processing  is  one  of  continuous 
evolution  of  techniques,  and  it  is  unrealistic  to  assume  that  steady -state  has  been  reached. 

The  preparation  of  reports  in  this  series  is  under  the  Editorship  of  Charles  R. 
Berndtson  of  Lincoln,  and  of  D.  Nessman  and  R.  French  of  Philco-Ford  Corporation. 
Inquiries,  suggestions,  corrections,  criticisms,  and  requests  for  additional  copies 
should  be  directed  to  C.  R.  Berndtson. 

The  principal  contributor  to  this  report  was  G.  L.  Shapiro  (Philco-Ford). 

Due  to  the  intricate,  evolutionary  manner  in  which  the  programs  came  into  being,  f  u: 
editors  regret  that  it  is  in  general  impossible  to  give  due  credit  to  all  —  mathematicians 
or  radar  analysts  or  programmers  --  who  contributed  to  the  definition  and  writing  of 

the  programs. 
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COMMON  SYMBOLS  AND  ABBREVIATIONS 


(The  units  given  for  certain  quantities  are  the  units  commonly  used  for  those 
quantities,  unless  otherwise  noted. ) 


ADT  ALCOR  Data  Tape 

ALCOR  ARPA -Lincoln  C-band  Observables  Radar 

ALTAIR  ARPA  Long-Range  Tracking  and  Instrumentation  Radar 

Alt  Altitude  (km) 

APS  Average  Pulse  Shape 

ARS  ALTAIR  Recording  System 

Avg  Average,  Averaging 

Az  Azimuth  (deg) 


c  Speed  of  Light 

CADJ  Adjusted  Calibration  Constant  (db) 

C-band  ALCOR  frequency,  5664  MHz  (NB)  and  5667  MHz  (WB) 

DBLT  Wide  Band  Pulse  Doublet 

El  Elevation  (deg) 

EOF  End  of  File 

GMT  Greenwich  Mean  Time 

h  Hours 

Hz  Hertz 


IF  Intermediate  Frequency 

in  Inches 


LC  Left  Circular  Polarization 

Isb  Least  Significant  Bit 

min  Minutes 


NB 

NRTPOD 

POD 

Phase 

PRF 

PRI 

pps 

pts 


Narrow  Band 

Non -real  Time  Precision  Orbit  Determination  Program 

Project  PRESS  Operation  and  Data  Summary  Report 
Presented  in  deg 

Pulse  Repetition  Frequency  (pps) 

Pulse  Repetition  Interval  (s) 

Pulses  per  second 
Points 


R 

• 

R 

rad 

RC 

RCS 

RF 


Range  (km) 

Range  Rate  (km/s) 

Radians 

Right  Circular  Polarization 
Radar  Cross  Section  (dbsm) 
Radio  Frequency 


s 


SDBLT 

S/N 

T 

TAL 

UHF 

V 


V 

w 

VHF 

WB 

WBS 

WTR 


e 

\ 

* 


Seconds 

Standard  Deviation  of  Wake  Velocity 

Wide  Band  Slaved  Pulse  Doublet 
Signal -to-noise  Ratio 

Time 

Time  After  Launch  (s) 

ALTAIR  Frequency;  415  MHz 
Velocity 

Doppler  Velocity 

Mean  Wake  Velocity 

ALTAIR  Frequency;  455.5  MHz 

Wide  Band 

Wide  Band  Slaved 

Western  Test  Range 

Total  Off-axis  Angle  (deg) 
Wavelength 

Denotes  Multiplication 


vii 


FLOW  DIAGRAM  SYMBOLS 


PROCESS,  ANNOTATION 


DECISION 


TERMINATOR 


SUBP'V’TINE:  where  NAME  le  the  entry 
call  Into  the  eubroutlne 


CONNECTOR:  where  P  specifies  a  page  In  the 

flow  diagram,  and  L  designates 
a  statement  number  In  the  program 
listing  or  a  refemce  point  In  the 
flow  diagram 

CONNECTOR:  where  X  Implies  a  continuation 
of  the  diagram  to  the  next  page 


INPUT/OUTPUT  OPERATION 


MAGNETIC  TAPE 


PUNCHED  CARD 


DISK 


ALCPOD 


I.  PURPOSE  AND  UTILIZATION 

A.  Source  of  Data 


B. 


ALCOR 
Data  Input 

ALCOR  Data  Tape  (ADT) 

Ascription 

n  An  ALCPOD  is  designed  to  produce  punched  card  metric  data  on  hard 
body^targets^in  a  format  suitable  for  input  toJ^RTPO&.  The  data  are  ultimately  used 
to  obtain  a  £  profile.  ALCPOD  is  norrpaliy  run  every  0.1s  without  smoothing. 

Output 

1. 


A 

D. 


2. 


/ 


A  listing  of  all  computed  quantities. 

Punched  cards  containing  corrected  R,  Az,  and  El  in  a 
format  suitable  for  input  to  NRTPOD.  (Optional:  smoothed  R. ) 


/ 


II. 


DESCRIPTION 


ALCPOD  produces  punched  metric  data  for  input  to  NRTPOD.  The  program 
processes  either^NBor  WB  tapes,  and  uses  only  primary  pulses. ^The  time  between 


cards  is  determined  by  NSKIP,  the  number  of  primary  pulses 
Therefore,  the  sampling  rate  changes  as  the  PRF  changes. 


en  output  samples. 


GMT  (h,  min,  s)  is  derived  from  one  of  three  sources: 

a.  The  transmitted  pulse  time  (accurate  to  10  ps)  from  tape. 

b.  Calculated  time  using  the  PRI  and  the  received  time  of  the 
initial  pulse  (accurate  to  1  ms). 

c.  The  received  pulse  time  (accurate  to  1  ms)  from  tape. 

R,  Az,  and  El  are  corrected  as  follows: 

R  =  IRANGE  +  TRBIAS  +  TTCOR  +  RRCOR  -  RCORF 

Az  =  IAZ  +  A  Z  BIAS 

El  »  IEL  +  ELBIAS  -  ECORF 

where 

IRANGE  is  uncorrected  R 

TRBIAS  is  range  bias 

TTCOR  (transmit  time  correction)  =  RR/c 

RRCOR  is  range  doppler  coupling  correction 

RCORF  is  tropospheric  refraction  correction 

IAZ  is  Az  encoder  angle 

AZBIAS  is  Az  bias  (Calibration  Record  Word  602) 

IEL  is  El  encoder  angle 

ELBIAS  is  El  bias  (Calibration  Record  Word  603) 

ECORF  is  tropospheric  refraction  correction 
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R,  if  requested,  is  computed: 

N 

.  12  ^  [nR  (t  +  nAt)  -  nR  (t  -  nAt)] 

R(t)  = - 52 - - - - 

(At)  (2  N+  1)  (  (2N+1)  •  1] 

where 

R(t)  =  Range  rate  at  time  t 
N  =  No.  of  points  in  a  half  interval 
R(t  +  nAt)  =  Range  at  time  (t  +  nAt) 

At  =  time  between  consecutive  range  samples;  i.e.  l/(primary 
pulse  PRF) 

Before  processing,  the  main  program  checks  that  ITBAND  (tape)  =  IBAND 
(input).  This  determines  that  if  WB  data  is  requested,  WB  data  exists  on  ti^e  tape 
requested. 


III. 


OPERATION 


A.  Input 

Waveform 

Time  of  pulse  option 
R  option 

No.  of  pulses  in  smoothing  half  interval 
Skip  interval  (pulses) 

First  and  last  pulse  nos.  of  processing  intervals 
No.  of  processing  intervals 

A  sample  input  is  shown  in  Appendix  A. 

CARD  1  (515,  2X,  A3) 


(Col.) 

1-  5 

NBAND 

0  =  NB;  1  =  WB 

6-10 

ITLUSE# 

0  =  transmitted  pulse  time  accurate  to  10  ps  (0) 

1  =  calculated  pulse  time  accurate  to  1  ms 

2  =  received  pulse  time  accurate  to  1  ms 

11-15 

NVALS 

NINT# 

No.  of  processing  intervals 

16-20 

0  =  R  on  cards  (0) 

1  =  no  R  on  cards 

21-25 

NSMTFT 

u 

TITL 

No.  of  pulses  in  smoothing  half  interval  (5) 

28-30 

3  character  title  (ALC) 

CARD  2 

(6110) 

(Col.) 

1-10 

NSTART(l) 

First  pulse  no.  of  initial  processing  interval 

11-20 

NSTOP(l) 

Last  pulse  no.  of  initial  processing  interval 

21-30 

NSKIP(l) 

No.  of  pulses  to  skip  between  each  output  card 

31-40 

NSTART(2) 

First  pulse  no.  of  second  processing  interval 

If  left  blank,  program  s^s  to  indicated  value. 
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(Col.) 

41-50  NSTOP(2)  Last  pulse  no.  of  second  processing  interval 
51-60  NSKIP(2)  No.  of  pulses  to  skip  between  each  output  card 
Repeat  card  2  as  necessary. 


B.  Output 


LISTING 

Selected  input  parameters 
GMT  h,  min,  s,  and  ms 
R,  R,  Az,  El,  and  pulse  no. 


PUNCHED  CARDS 
Radar  identification  (A3) 


Year 

(3X,  12) 

Month 

(12) 

Day 

(12) 

h 

(12) 

min 

(12) 

s 

(12) 

fraction  of  second  (13,  3X)  or  (15,  IX) ^ 


Orbit  no.  (II) 

Az  (F8.3) 

El  (4X,  F8.3) 

R  (4X,  F12.4) 

R  (if  requested)  (5X,  F8.5) 


Sample  outputs  are  given  in  Appendix  B. 


IX)  if  ITLUSE  =  0;  (13,  3X)  if  ITLUSE  =  1  or  2 
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IV. 


PROGRAM  LIMITATIONS 


NVALS 
Length  of  run 


^  50  processing  intervals 

If  R  is  desired,  no.  of  pulses  :?  4000 
If  R  is  not  desired,  no  limit 


If  R  is  requested  (NINT  =  0),  a  processing  interval  should  not  contain  a  PRF 
change.  Program  processes  primary  pulses  only. 


V. 


PROGRAMMING 


A.  ALCPOD  (see  Appendices  C  and  D. ) 

ALCPOD  is  the  control  section  of  ALCPOD.  ALCPOD  reads  the  input 
cards,  calls  READJS  and  UNPACK,  and  punches  the  desired  metric  uata. 

B.  HEDADT  (see  Appendix  E. ) 

Subroutine  HEDADT  unpacks  the  ADT  header  record  which  contains 
bandwidth,  reel  no.,  WTR  no. ,  data  of  mission,  and  mission  designator.  The  call 
statement  is  HEDADT  [ISIG,  INBUF(l),  IEQM(l)  ] 

INPUT 

INBUF(l)  First  word  in  the  ADT  header  record^ 

OUTPUT 


IEQM(l) 

IZBAND 

(bandwidth:  1=W3,  0=NB) 

IEQM(2) 

ITREEL 

(reel  no.) 

IEQM(3) 

ITWTR 

(WTR  no.) 

IEQM(4) 

IMTH 

IEQM(5) 

IDAY 

(Dpte  of  test) 

IEQM(6) 

IYR 

IEQM(7-9) 

ITDESG 

(mission  designator) 

C.  UNPACK  (see  Appendix  F.)2 

Subroutine  UNPACK  unpacks  the  raw  data  from  the  ADT,  and  trans¬ 
lates  it  into  a  format  usable  by  the  IBM  360/67  computer. 


Not  used. 

i 

INBUF(2)  to  INBUF  (1803)  contain  the  remaining  yords  in  the  record. 
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D.  READJS2 

The  first  call  to  subroutine  READJS  opens  the  file  and  reads  the  ADT 
header  record.  The  second  call  to  READJS  reads  the  ADT  calibration  record  and 
stores  the  values  in  a  buffer  area.  ALCPOD  extracts  the  individual  calibration  values 
it  requires.  Each  subsequent  call  to  READJS  reads  an  ADT  data  record  consisting  of 
eight  ALCOR  pulses. 

E.  TIMDP  (see  Appendix  G.) 

TIMDP  converts  GMT  total  s  to  h,  min,  s,  and  ps. 

The  call  statement  is  TIMDP  (TIME,  IHR,  MIN,  ISEC,  IF  RAC). 

INPUT 

TIME  Time  of  pulse  transmission  (GMT  total  seconds  in 


double  precision) 

OUTPUT 

IHR 

Hours 

MIN 

Minutes 

ISEC 

Seconds 

IFRAC 

ps 

F.  SMOOTH  (see  Appendix  H.) 

Subroutine  SMOOTH  computes  R  using  the  original  R  and  T. 
The  call  statement  is  SMOOTH  (N,  L,  X,  NO,  ZH). 

INPUT 

N  No.  of  points  in  smoothing  half  interval 

L  No.  of  pulses  within  all  processing  intervals 

X  Array  of  Ranges 
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NO  Code  specifying  operation  desired.  Set  to  zero  to  obtain  R. 
ZH  Time  (s)  between  range  samples  (1/primary  pulse  PRF) 

STORED  IN  COMMON 

RR  Array  of  ft 

G.  DREFC  (see  Appendix  J. ) 

The  tropospheric  refraction  correction  subroutine,  DREFC,  is  based 
on  tropospheric  refraction  tables  in  PPP-36.  ^  A  modified  version  of  this  subroutine 
is  now  in  use.  DREFC  is  the  same  as  REFC  except  that  the  values  or  constants  ar^ 
defined  as  double  precision  words. 

The  call  statement  is  DREFC  (E,  R,  DEE,  DRR). 

E  =  Uncorrected  El  (must  be  between  0°  and  90°) 

R  =  Uncorrected  R  (km) 

DEE  =  El  tropospheric  correction 
DRR  =  R  tropospheric  correction  (km) 

The  corrected  values  to  be  computed  after  exiting  from  the  DREFC 

subroutine  are: 

El  =  E-DEE 
R(km)  =  R-DRR 


REFERENCES 
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appendix  b 

ALCPOD  OUTPUTS 


PUU-ALCOR  POLAR  *  LC  8AN0  *  NU 

ITLOSc  =  0 

RfcfcL 

NO.  « 

l 

TITLE  * 

ALC 

OATfc  * 

3/  3/7i 

START  STOP 

27191  27691 

SKIP  START 

9 

STOP 

SKIP 

START 

STOP 

SKIP 

START 

STUP 

S<  (P 

musE  =  o  :  ust  air  as  computeo  sv  double  prec.noru'  on  tape  utr-tti/21 

I T L USE  *  1  '  USE  GMT  AS  COMPUTED  8VPR0GRAM  AS  A  FUNCTION  OF  PRF  (TIME  OF  RECEFTIONI 

ITLUSE  *  2  :  USE  GMT  AS  COMPUTED  8Y  REGULAR  GMT  MOROS  (TIME  OF  RECEPTIONI 


MINT  =  0  (SMOJTHING  IS  TO  OCCURI 

NO.  OF  PO(NTS  ON  EACH  SIDE  OF  I  NT  EKPOL  AT  Ml  POINT  ■  5 


Kt  AL>y 

TU 

SMOOTH 

501  1 

PUINTS  WITH 

A  DELTA  T 1  ME ( SEC  I 

*  0.310 

THEN  PJNCrt  EV 

ERY 

I  ITLE 

DATE 

H  M 

S  MS 

AZIUEGI 

ELI  JEGI 

RANGE ( K  M 1 

ROOTIKM/SI 

ALC 

71 

3 

3 

544 

7. 30910 

0 

61.860 

19. 796 

569.7701 

0.00017 

27191 

ALC 

71 

S 

3 

544 

7.40910 

0 

ot.571 

19.  760 

569. 1C02 

-6.69679 

27201 

ALC 

7l 

3 

3 

544 

7.50910 

0 

61.576 

19.762 

563.5317 

-6.68650 

27211 

ALC 

71 

3 

3 

544 

7.60910 

0 

01.  587 

19. 768 

562.7609 

-6.71191 

27221 

ALC 

7i 

3 

3 

544 

7.70909 

0 

61.595 

19. 7o6 

562.0908 

-6.69662 

2  72  31 

ALC 

71 

3 

3 

544 

7. 80909 

0 

61.605 

19.771 

561.9219 

-6.69025 

2  7241 

ALC 

71 

3 

3 

544 

7. 90909 

0 

61.  606 

19.773 

560.  7525 

-6.69563 

2726i 

ALC 

7 

3 

3 

344 

b.  909 

0 

61.609 

19.782 

560.0829 

-6.69636 

2  7261 

ALC 

7. 

3 

3 

344 

8. 10°0d 

0 

oi.612 

19.785 

559.5137 

-6.69297 

2  72  71 

ALC 

71 

3 

3 

544 

8.20908 

0 

61.615 

19.790 

558. 7538 

-6.69912 

27281 

ALC 

71 

3 

3 

544 

8.3  0908 

0 

6 1  .62  0 

19.793 

558.0752 

-6.69620 

2  72  91 

ALC 

71 

3 

3 

544 

8. 4090b 

0 

61.623 

19.795 

557.5055 

-6.69766 

27301 

ALC 

7i 

3 

3 

544 

8.50907 

0 

61.626 

19.801 

556. 7397 

-6.69755 

273H 

ALC 

71 

3 

3 

544 

8.60907 

0 

01.628 

19.809 

556.0696 

-6.70117 

27321 

aLC 

71 

3 

3 

544 

8. 70907 

0 

61.631 

19.  806 

655.3997 

-6.69925 

27331 

ALC 

7l 

3 

3 

5*4 

8.80907 

0 

61.639 

19.815 

555.7248 

-6.69957 

27351 

ALC 

7i 

3 

3 

344 

H. 90907 

0 

61.637 

19.817 

554.0548 

-6.69996 

27351 

AL  C 

71 

3 

3 

544 

9.  906 

0 

6 1.  692 

19.826 

553.3856 

-6.70181 

27361 

ALC 

71 

3 

3 

544 

9 •  109  3  o 

0 

61.698 

19.828 

552.7157 

-6.69980 

2  73  71 

ALC 

71 

3 

3 

544 

9.  2C906 

0 

61.698 

19.837 

552.0553 

-6.59862 

27381 

ALC 

71 

3 

3 

544 

9. 30906 

0 

61 • 650 

19. 839 

551.  3757 

-6.70182 

27391 

ALC 

71 

3 

3 

544 

9.40905 

0 

61.659 

19.895 

550.7056 

-6.70105 

2  7501 

AL  C 

71 

3 

3 

544 

9. 50905 

0 

61.661 

19.  898 

550.0355 

-6.70179 

27511 

ALC 

71 

3 

3 

544 

9.60905 

0 

61 • 669 

19. 850 

559.3651 

-6.70333 

2  7521 

ALC 

71 

3 

3 

544 

9.  70903 

0 

61.670 

19.856 

558.6939 

-6.70253 

27531 

ALC 

71 

3 

3 

544 

9. 80904 

0 

61.672 

19.  8  62 

558.0236 

-6.70257 

27551 

ALC 

71 

3 

3 

344 

9.90904 

0 

61.678 

19.867 

557.3531 

-6.70485 

27451 

ALC 

7 1 

3 

3 

5441 C.  904 

0 

61.063 

19.873 

556.6328 

-6.70332 

27461 

ALC 

71 

3 

3 

54410. 10904 

0 

61.692 

19. 881 

556.0125 

-6.70382 

27471 

aLC 

71 

3 

3 

54410.20904 

0 

61.699 

19.889 

555.3515 

-6.70907 

27481 

ALC 

71 

3 

3 

54410. 30903 

0 

ol. 697 

19.889 

554.6710 

-6.70552 

27491 

ALC 

71 

3 

3 

54410. *0903 

0* 

61.703 

19.892 

544.0005 

-6.7  0582 

2  7501 

AL  c 

71 

3 

3 

34410.50903 

0 

61. 706 

19.897 

543.3300 

-6.70529 

2  7511 

ALC 

71 

3 

3 

5441C. 60903 

0 

61.  715 

19.  903 

542.6567 

-6.71180 

27521 

«L  C 

71 

3 

3 

34410.70902 

0 

61.716 

19.906 

541.9881 

-6.70636 

27531 

ALC 

71 

3 

3 

5*410.  80902 

0 

61.719 

19.911 

541.3177 

-6.70563 

7541 

ALC 

71 

3 

3 

54410.90902 

0 

61.  725 

19.919 

550.6570 

-6.70725 

27551 

mLC 

?4 

j 

3 

5*411.  902 

0 

61.730 

19.922 

539.9763 

-6.  10721 

2  7561 

ALC 

i  1 

3 

3 

54411.  10901 

0 

61. 735 

19.925 

539.3356 

-6.70753 

27571 

ALC 

71 

3 

3 

54411.20901 

0 

61.738 

19.930 

538.6358 

-6.70739 

27581 

«LL 

71 

3 

3 

54411.30901 

0 

61.791 

19.936 

537.9651 

-6.70793 

2  7591 

ALC 

71 

3 

3 

5*411.40901 

0 

61. 796 

19.939 

537.2932 

-6.70872 

27601 

ALC 

71 

3 

3 

54411 .50900 

0 

61.752 

19.959 

536. 6225 

-6.70855 

27611 

12 


ALC  71  3  3  544  7.40910  0  61.571  19.760  564.1002  -6.69879 

00000000000000000000  00  0  0000000000000000000  0000000000000  0000000000000000000 

1  2  3  J  S  G  I  8  9  10  H  12  13  14  *5  16  12  18  19  20  21  22  23  28  25  26  27  28  29  3031  32  33  34  35  36  32  38  39  40  41  <2  43  84  85  <682  88  88  50  51  52  53  58  SS  56  57  58  59  60  Dl  62  63  68  05  66  62  68  69  20  21  22  23  24  25  26  11 28  29  80 

min  11111111111111  min  m  mm  mmimimm  iimimmimmii 

2  2222222  222  22222 122222222222222222222222 2 22  2  222222222222222 22  222121222211222122 

0  333333  3  333333  33333333333  33333333333  33333333333333  33333333333  3333333333 

^4444444444444  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  H  4  4  4  4  4  4  4  4  4  1  4  4  4  i  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4 

5555535555555  55555555555555555  5555555555555505555555  5555555555555555555555555 
666fi6C6C66606666G66666666666  6 6 6 6 6 6 6 6 6 6 6 6 G 6 6  6666666666  666666666666  6  666666666 
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APPENDIX  C 


ALCPOD  PROGRAM  LISTING 


DOUBLE  PRECISION  RDOT ,  OTRB, RANGE • T TCOR  t DR RUSE • RRCOR  *  LI »XAZ, DAZBt 

1ZL ,XEL  »  DELB ,RNGF , EL VF , R ADEL  »ECGRF » RCORF , XCPT IM 
RE AL*8  KANG.RR 

DIMENSION  IYEARUOOO) .MONTH  I  4000  )  ,  I  DAY  1 4000) . 1  o  ?°  c  .*  unom 

1  M INI 4C00 ) *  I SE C 1(4000 1* I  SEC  2 1 4000 ) t  AZ  ( 4000 ), EL ( 4000 ) , 

2  KANG(40C0)«RK(4CCC)« ISPRII4000) 

DIMENSION  NSTA'  T I  50  I • NSTOP I 50 > » NSK I P ( 5C 1 
DIMENSION  XATBK 1 28  I , Oe I  AS  I  8 »  IPRSI8I , lADD(fl) 

DIMENSION  XNBUF  I  18C3I  »  P  IFA(  ib  I  tOIF  AH6  l  »  XKR-SI5) 

DIMENSION  IECMI9) t IIDESGI  3) 


IXDPTIM* ICAT (68?) 


C0UIVALENCE(XNPUF(  l»,  INBUFIl)  »  ITDCC1.  ,tcnym  itutri 

EQUIVALENCE  (  I FQM (l» , IZBAND » , ( IEQM I  2 ) » I TREED , IIEQM 1 3 ) , 1 TW TR ) . 
? ( I EGM ( 4  )  ,  I  TMNTH) »  UECM  (  5  )  , 1 TDAY ) ,  (  I  EGM  6  ) » I TYE AR) • 

3  (  I  EGM  (  7)  ,  I  T  D  E  SCI  {  1  )  ) 

ZLCZ’LC  '/*ZRC/*RC  '/,ZWB/’WB  »/»ZNB/*NB  •/  ,NTCT/0/ 

IFRST3/0/, IFRST4/0/,lNTAV/l/» IFRST2/0/ 

ALC/ ' ALC* /, ICN/I/t IZERO/O/i IFRSTl/O/tIMSAVE/O/ ,BLNKK/» 
IPRS/20C, 160,100,80,50,40,25,20/ 

I  ADD / 1 0 1 1 7 , 10* 13,20,25,40,50/ 


DATA 

DATA 

DATA 

DATA 

DATA 


I  PUL AR 
I  POLAR 
NBAND 
NUAND 
ITLUSE 
l  TL  JSE 
ITLUSE 
NEWA  = 
NEWA  = 
NINT  = 
NINT  = 


=  0 
=■  l 
=  2 
C 

1 
0 
I 


*  o  LEFT  CIRCULAR  DATA  REQUESTED 
=  1  RIGHT  CIRCULAR  DATA  REQUESTED 
=  0  NARROW  BAND  DATA  REQUESTED 
:  I  WIDE  BAND  DATA  REQUESTED 

LSE  GMT  AS  COMPUTED  BY  DOUB . PREC.TP.WD. ( ( Tfc-TT ) 12  ) 
USE  GMT  AS  COMPUTED  BY  PRF  IN  PROGRAM  (TR) 

USE  GMT  AS  COMPUTED  BY  REGULAR  GMT  WCRCS  ITR) 
MISSION  FLOWN  BEFORE  15  CCT  70  (OLD  ATTN-) 

MISSION  FLOWN  AFTER  li  CCT  70  INEW  ATTN.) 

SMOOTHING  IS  DONE 
NO  SMOOTHING 


RE  AC ( 5 •  I )  NBAND » ITLUSE  tNVALSt NINT tNSMOOl.T I TL, 

2 ( NST  AR  T (  I  ) , NSTCP ( l ) ,NSK I P  (  I  ) , l*l,NVALS> 

1  FORMAT (515,  2X,A3/(6I10)) 

IF (NSMCOI.LE.O)NSMCCI=5 
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NSMOG  =  N  SMOG  1 
MSVE=NSK I P ( 1 ) 

IF<TITL.NE.BLNKK»ALC»TJTL 
NSW= I TLUSE 
C 

IE0F*U 
I  ERR  =  0 

CALL  REACJSI INPUFV IEOF, IERR) 

IS IG“  1 

CALL  HEUACT  ( I S IG , I NBUF ( I ) , I ECM (1) ) 

I T  BAND3 1 1 BAND 
NEWA=0 

ir(ITYEAR.GT.7DGC  TO  282 
IF(ITYEAR.LT.?DGG  TO  283 
IF<  ITMNTH.GT.iDGO  TO  282 
IFdTMNTH.LT.lDGO  TO  283 
I  F  (  ITDAY.LT.  1 I GC  TC  283 

282  NEWA= 1 

283  CCNT  INGE 
I  FRR“0 

CALL  REACJSI INPUF, 1ECF, IERR) 

C 

C  STORE  THE  CFSIRED  CALIBRATION  VALUES 

C 

N=0 

DU  20  K*256»  387 
N=N+  I 

20  XATBL(N)*XNBUF(K) 

C 

N*0 

DO  22  K=5I2i 527 
NXN+  I 

22  P  1  F A ( N I =XNBUF I K ) 

N=0 

DO  7 3  K*528,5<t1 
N*N+  I 

23  O1FA<N)=XN0UFIK) 

C 

PFSA*XNBUF<592> 

PSSA=XNBUF<  5931 
DFSA=XNBUF ( 59A I 
0SSA=XNBUF<595) 

C 

APIAS=XNBUF(60?  ) 

EblAS«XNBUFI60D 

DEGCON= ( 18C.A.0A79369I/31A1.59 

azbias=degcon*abias 

ELBIAS*OEGCCN*FBIAS 

C 

N  =  0 

DO  25  K  =  604,6U 
N*N*  I 

QB  I  AS ( N ) =XNBUF ( K  ) 

25  RBIAS(NI=GBIAS(N) 

C 

PWRCN=XNBUF(62D 
PWRSN  =  XNBUF(621  ) 
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—if  riniMMUiiiT  ii  Hi ii  u  —  - . 


PWKCW=XNBUF(62?1 
PWRSW  =  XNBUF(62',1 

C 

N*0 

DO  27  K  =624 , 62P 
N  =  N»  1 

27  XKRCS(Nt=XNBUF(K) 

C 

PSSL=XNBUF(629» 

GSSL=XNBUF!6301 

C 

JCL'N  =  -  1 
INDEX  =  (J 
I T  ST* 1 
I  TDEf  =  1 
IPOLA«*0 
I TCNT=0 
JJ=0 
IPULS=0 

c 

DC)  120  I J*  1  ,N VALS 
NBEG=NSTART ( IJI 

C 

NNSET-NSKIP! I J I ♦ 1 
IF(NINT.EC.0)NNSET*1 

c 

3  JCCM=JCCN+1 

IF(  JCON.EC.9.0R. JCCN.EC.C1GC  TO  97 
INDEX=<  JCCN-1»*9C0 
GO  TO  99 

97  JCCN=1 

I  N  D  E  X  =  0 

98  I  EOF  =  0 
I ERR  =  0 

CALL  REACJS!  INPUF, lEOf ,IERR1 
IF ( IERR.EG. 1 1G0  TO  103 

99  CA^L  UNPACK 

IF!  IC0CE.EC.3.C]R.ICCCE.EC.7.0R.IC0CE.EC.21GC  TC  620 
I F ( ITLUSE.NE.01GC  TO  ICO 
CALL  TIMCP(XCPTlM,IHRS,iMIN,lSEC, I  MSEC  1 
100  CONTINUE 

IF( IFR5T2.E0.  1IGC  TO  92 
ZBAN=ZNB 

IF ( I TBANC.EQ. l»ZBAN*ZWB 
ZPOL=ZLC 

IF(  I  PULAR.ECI.il ZPOL*ZRC 
RRUSE*-. 00943 

I F ( I TB AND. EQ.11RRUSE*-. 000115 

C 

WRITE(6,200»ZPnL,ZBAN,ITREEL,TITL, IIECP(I), 1*4,61 
200  FORMAT (' 1P00-ALCCR  POLAR  ■  • , A2 , 4X, • BANC  =  * , A2,4X, ' REEL  NO. 
1,15,'  TITLE  =  ' , A4 ,  '  DATE  *  • , 1 2, * /• , 1 2 , • / • , 12 1 
WR  I TE  I  6 , 208  1  ITLUSE 
208  FORMAT!  •  ITLUSE  *  ',131 

WRITE (6 ,2 12  > (NSTART! I ),NSTOP(  I),NSKIP( I  1 , I  * l , NVALS 1 
212  FORMAT ( • 0  START  STOP  SKI P* , 10X  START  STOP  SKIP'.ICX, 

1  'START  STOP  SKIP' , 10X, 'START  STCP  SKIP'/ 
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2  <4<2X,  1 5  i  2X  i  I  5 , 2X  ,  I5.PXJ  )  ) 

WR ITE ( 6 , 20 1 ) 

201  FORMAT! 'OITLUSf  =  0  0  USE  GMT  AS  COMPUTED  BV  DOUBLE  PREC. WORDS  ON 
ITAPE  ( ( TR-TT )/? »  •/•  ITLUSE  «  1  0  USE  GMT  AS  COMPUTED  BY 

2PR0GRAM  AS  A  FUNCTION  OF  PR?  (TIME  CF  RECEP T I  ON )  • / •  ITLUSE  «  2  0  U 
3SE  GMT  AS  COMPETED  BY  REGULAR  GKT  WORDS  (TIME  OF  RECEPTION)  •/) 

IF ( N INT.NE.O )GH  TO  2  IB 
WRITE(6,214)  NSMOO 

2l<i  FORMAT! 'CNINT  =  C  (SM0C1HING  IS  TO  OCCUR)1/  IX, 'NO.  CF  POINTS  CN  E 
IACH  SICE  CF  INTERPOLATED  POINT  ■  *,14) 

218  CONTINUE 

IF (NBAND.NE. ITPANQ)GO  TO  695 
l  FRST2= l 
92  CONTINUE 

IF!  ITLUSE. NE. I JGO  TU  619 
IFINUMPRI.LT. NSTARTII JIIGO  TO  6i6 
IF!  IFRSTl.EQ.DGC  TO  617 
I  ACMS=0 
DO  612  K- 1 ,8 

612  IF!  IPRS!K).EC.IPRF)IADMS*IADD(K) 

IF!  IAOMS.GT.CJGO  TO  602 

599  WRITE! 6,601 )NUMPRI , IPRF 

601  FORMAT! 'O' , 12X, 'UNKNOWN  PRF • , 3X ,  •  PR  I .  NO  -  •  ,  I  10 , 5 X  ,  • PRF  *  •  , I  5/ ) 
IF!NINT.EC.O)Gr  T  !]  121 

GO  TO  125 

602  IHNEXT.IHRS 
I NNEX  T  = IMIN 
I  SNEXT  = I  SEC 

584  I MNEXT  = I  MSEC  +IACPS 

IF!  IMNEXT.LT. IOOOGO  TC  58‘ 

IMNEXT=IMNEXT-10C0 
l SNEX  T3  I  SNE  XT* 1 

IF  (  ISNEXT.LT. 6OG0  TC  585 
ISNEXT=ISNEXT-60 
I NNEXT* I NNEXT  + I 

in  INNEXT.LT. 60  JGO  TO  585 
INNEXT*INNEXT-60 
I HNEXT  =  IHNEX  T ♦ 1 

585  I F  RST l *  1 
GO  TO  616 

617  I  A  CMS=0 

DO  613  K=l, 8 

IF ! IPRSIK ).EC. TPRF  )  I ADMS* l ADD ! K  ) 

613  CONTINUE 

IT!  I  ADMS. Ltf. OHO  TO  599 

586  IHRSMHNEXT 
IM IN= INNE  XT 
I SEC= I SNEXT 

I M  SEC= I MNE  X  T 

I MNEX  T* I  MSEC  +IACMS 

IF!  IMNEXT.LT.  ICCOGC  TO  616 
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IMNEXTMHNEXT-10C0 

ISNEXT= ISNEXT*! 

IFIISNEXT.lt. 6H)GC  TO  616 
ISNEXT=ISNEXT-60 
INNEXTMNNEXT  +  ! 

IF(  INNEXT.LT. 6CIG0  TO  616 
INNEXT= INNEXT-60 
IHNEXT= IHNEXT* 1 

L6  CCNTINLE 
19  CONTINUE 

10  I F ( NUMPR  I.LT.NSTART ( I J ) )GO  TO  3 
I F ( NUMPR  I.GT . NS  T  AR  T (  IJ  )  )GC  TO  628 

GU(T00627EC‘3,PR*ICCDE,EC“7*OR,ICODE*EC*2>GC  T0  600 
0  NSTART (  IJ)=NSTART (  l  J  )  ♦  1 
WRITE(  6,6314)  I  J,  NS  TARTU  J» 

'  e°h;"t5,onstakt‘’*>3,»)  HAS  BEEN  CHANGED  TO  SIIO) 

7  I  PUL  S  =  0 
GO  TO  10 

8  lpi<Is  =  5pULtS;^CR*,COOE*EC-7-°R-ICODE-EC-2,GC  T0  118 

IF (  IPULS.NE.NNSET )G0  TO  118 
IPULS=0 

D  CONTINUE 

IF  (  IFRST4.NE.01GC  TO  11 
DT  IME*  (  1. /FLOAT  (  IPRE  J  ) 

1  ”?S^R0GT??J^!DCT,/,8l92*OD"OO,,*l4'-9e9625D^00 

OTRH=TRB IAS 

?muSEA:;,T5c^ 

IDELTM=(RANGF/?99776.D*00)*1.0D*06 

!I!^ui:c^LGATn''SEC*IDELTW,/lo*OD+oof.5r+oo 
IF( IMSEC.LT. lOOOCOJGO  TO  641 
IMS EC-  IMSEC-IOOOCO 
ISF.C=ISCC*1 

IF( ISEC.LT. 60100  TO  641 
ISEC=ISEC-60 

imin=imin*i 

CONTINUE 
GO  TO  664 

II5r?=<RANGE/299776*D*co,*,ROOT/i°oc.oc+oo) 

'lANGfc  =  RANGE*TTrOR 

CONTINUE 

DRRUSE=RRUSE 

RRCOR=ORRUSE*ROOT 

9ANGE=RANGE*RRC0R/ 1C00. D*00 

0AZ8=AZBIAS 
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XAZ=XAZ*OAZB 

Z |_  —  I  OFLOAT  (  I  EL)  *2. 0*00*3141. 59265350*00/12.  OD+OO**  17) 

XEL=ZL*. 057295*0*00 
OELB  =  ELB I  AS 
XEL=XEL*CELB 

CALL  DREFC  ( XEl  , RANGE , ECCRF ,RCQRF ) 

RNGF=RANCE-RCOPF 
ELVF=XEL-ECQRF 
RADEL=LLVF*. 0174530*00 
RANGE=RNGF 

IFININT.EC.O)  IPN*0 

IF! ION*  EG • l)WRITE(6,647S 

647  FORMAT! ’OTITLE  DATE  H  M  S  MS  AZ(OCG)  EL(CEG)  RANG 

IF  i  KM »  RCCKKM/S)  »/) 

I  0N=0 

IF (NINT.NE.O)GO  TO  681 

NT OT  =  NTGT*  l 

IYEAR(NTOT)=lTYEAK 

MONTH! NTOT ) » I TMN  TH 

I  DAY ( NTOT )  =  I TOAY 

IHOUR!  NTOT  I  *  I  Ht»  S 

MIN(NTCT)*IMIN 

I SEC1 ( NTOT  1= I SPC 

I$EC2(NTCT)=IMSEC 

AZINTOTMXAZ 

EL ( NTO  T ) =EL  VF 

RANG! NTOT )=RNGT 

ispruntoti^numpri 

IF !NTOT.EQ.4COO)GO  TO  121 

GO  TO  118 

681  IF ! ITLGSE.NE.O)GC  TO  117 

WRITE! 6  ,645) ALC » I TYE AR , I TMNTH , I  TO AY , I  HRS, I  MIN, I  SEC , I  MSEC , I  ZERO, 

IXAZ,ELVF ,RNGF,*ZCT,NGMPRI 

645  FORMAT ! IX , A  3 , 3X , 6  I  2 , » . » , I  5, IX , I l , F 8. 3 , 4X ,F8. 3 , 4X ,F 1 2.4 , 5X , F 8. 5  , 
1110) 

WR ITE! 7  ,644) ALC, I TYE AR , I TMNTH , I TOAY , IHR S, I  MI N, I  SEC , I  MSEC , I  ZERO, 

IXAZ ,ELVF  ,RNGF 

644  FORMAT! A3,3X,6T2, ’ . * , I  5 , IX , I l , F 8. 3 , 4X , F 8. 3 , 4X ,F 1 2. 4 , 5X , F 8. 5 ) 

GO  TO  118 

117  WRITE! 6  ,650 ALC, I T YE AR , I TMNTH , I TD A Y , I  HRS , I  MI N, I  SEC , I  MSEC , I ZE RO , 

IXAZ»ELVF»RNGF»RZCT,NGMPRI 

65C  FORMAT! IX , A 3, 3X , 6 1 2 , » . *,I3,1X,1 1 ,F8. 3 , AX ,F 8 . 3 , 4X ,F 1 2. 4 ,5X , F8. 5 , 
l  I  10) 

WRITE! 7  ,649) ALC, I TYE AR , I TMNTH , I TOAY , IHR S , I M I N, I SE C , I  MSEC , I Z E RO , 

IX AZ  *  EL VF , RNGF 

FORMAT (A3,3X«612, » . » » I  3 , 3X , 1 1 ,F8. 3 ,4X , F8. 3 , 4X , F 12. 4 , 5X ,F 8 . 5 ) 

IF INUMPRI .LT.NSTCP! I J ) )G0  TO  3 
I PULS=0 
IFRST  1  =  0 
I FRST  3  =  0 
CCNT INGE 

IMNINT.NE.OGP  TO  125 
NUSE=MSVE*l 

WRITE! 6, 128) NTOT, DTIME.NUSE 


649 

C 

116 


119 

120 
C 

121 
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128  FORMAT ( • ORE ACY  TC  SMOOTH’, 16,'  POINTS  WITH  A  DELTA  TIMEISEC)  *  ’ 
1 , F6. 3, ’  THEN  PUNCH  EVERY  ’,16) 

WR I TE { t,647) 

NCARD$=NTOT 

CALL  SMOCTH(NSMOO,NCARDS,RANG,0»DTIME  ) 

I TESS*- 1 

DO  6  1=1 .NCARDS 

I T  ESS* I TESS*  1 

IK  ITLLSE.EQ.O )GC  TC  7C4 

IF  (MOD!  ITESS,NtlSE  »  .NE.C  )GO  TO  6 

WRITE! 6, 650) ALC, I  YEAR! I ) .MONTH! I ), I  DA Y(I) , I  HOUR! I ) , M I N (  I  ) , I  SEC  1 (I) 
1, I SCC2 ( I ) , IZERP.AZ! I ),EL( I ) .RANG! I ) ,RR II) ,  I SFR I  U  ) 

WRITE ! 7,649) ALT, I  YEAR! I ) .MONTH! I ) , I  DAY  1 1 ) , I  HOUR! I ) , MI N ( I ) ,  I  SEC  1 ( I i 
l, ISEC2! I ), IZERP.AZ! I ) ,EL( I ) .RANG! I ) ,RR( I ) 

GO  TO  6 
704  CONTINUE 

I  F  (  MOD  (  ITESS.NtJSEI.NE.O  )GC  TO  6 

WRITE (6,645) ALC, I  YEAR! I ) .MONTH! I ) , I  DAY ( I ) , I  HOUR! I  )  .MINI  I ) -  'SECl! I » 
1, ISEC2! I ), IZERr.AZ! I ) ,EL( I ) .RANG! I ) ,RR 1 1 ) , I SPR I (I) 

WR ITE! 7,644) ALC, I  YEAR ( I ) .MONTH! I ) , I  DAY! I ) , I  HOUR  I  I ),MINII),ISEC1(I) 
1, ISEC2! I ), IZERP.AZ! I ) .ELI  I ) , RANG ( I ) ,RR ( I ) 

6  CONTINUE 

GO  TO  125 

103  WRITE! 6, 107 ) NUMPR  I 

107  FORMAT! ’OPARITY  ERROR  ON  READ  AFTER  PRI  *  ’  .  I  1C ) 

GO  TO  99 

680  WRITE(6, 1091NUWPR1 

109  FORMAT!’  END  OT  FILE  REACHED  LAST  NUMPRI  VALUE  *  ’,110) 

GO  TU  125 

6' 5  ^RITE(6,114)NBAND,ITBAND 

114  FORMAT!'  INPUT  BAND*  ’I1C,’  BAND  ON  TAPE  =  *110) 

125  RETURN 
END 
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APPENDIX  D 

ALCPOD  FLOW  DIAGRAM 


READ 

INPUT 

CARDS 


INITIALIZE 

VALUES 


READIS 

READ 

HEADER 

RECORD 


- * - ^ 

HEDADT  \ 

UNPACK  \ 
HEADER  / 
RECORD/ 


I 


/REA  PIS 

READ 

\CALIBRATIOWy 
\  Record 


STORE 

CALIBRATION 

RECORD 
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’  'T’-sv” 


APPENDIX  D-2 


.J¥.  '  fV- 


v*r 


APPENDIX  D-3 
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if--- 


APPENDIX  D-4 


APPENDIX  D '5 


* 


XISiG 

XICAL 

XIEQU 

BASE 

HECACT 


RETURN 


TEMPI 

TEMP2 

SAVEA 


APPENDIX  E 

SUBROUTINE  HEDADT  PROGRAM  T.TSTTMr. 


CALL  HECACT  ( ISIG, IN6UE, I  ECU  > 

START  ,S,G  "  1  UNP*CK  ™E  20  k0RC  *DT  HEADER 

ENTRY  HEDADT 

SPACE 


EQU 

A 

ECU 

5 

ECU 

6 

EQU 

12 

SPACE 

SAVE 

(14, 17), T, A 

BALR 

12,  C 

USING  * • BASE 

ST 

13, SAVEAAA 

LA 

7,  SAVPA 

ST 

7,8(0,13) 

LR 

13,7 

SPACE 

LM 

XI  SIG.XIECU, 

0(1) 

SPACE 

L 

e.oixrcAL) 

ST 

8, TEMPI 

ST 

8, TEMP2 

SRL 

8,31 

ST 

B  « 0 ( X ! ECU ) 

MflANC 

L 

8, TEMPI 

SLL 

6,1 

SRL 

8,25 

ST 

a.A(XIECU) 

mreel 

SPACE 

L 

8 , A ( X ICAL ) 

ST 

8, TEMPI 

ST 

8, TEMP2 

SRL 

8,16 

ST 

8 , 8 ( X  T  ECU ) 

MfcTR 

L 

8, TEMPI 

SLL 

8,  16 

SRL 

8 , 2A 

ST 

8,  12(X|EQU) 

MMNTH 

L 

8,  TEMP2 

SLL 

8,24 

SRL 

8,2* 

ST 

8, 161  XIEQU) 

MCA  Y 

SPACE 

SR 

8,8 

1C 

8 , 8 ( X  I C AL  ) 

ST 

8, 2C ( X | ECU ) 

MYEAR 

MVC 

24(9,X|LCU),9(XICAL) 

SPACE 

L 

13, SAVE  A  +  A 

RETURN 

«  1 A ,  1  2  )  f  T 

CNCP 

0,  A 

CC 

F  ’  0  ’ 

CC 

E’O* 

DC 

18 A ( A  ) 

END 
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APPENDIX  F 

SUBROUTINE  UNPACK  PROGRAM  LISTING 


UNPACK 


START 


CUBUF 
BASA 
ST  ART  1 


C  SECT 


ENTRY  UNPACK 
SAVEL 
DROP  15 
CNOP  0,4 

BALR  2,0 
USING  START, 2, 3 
L  3, BASA 

L  4 ,  DUBIIF 

L  5 ,  DUB<  !F 

L  6 , OUBUF 

A  5 ,  *F  ’ 4  096  * 

A  6 , *F  *  P 192  * 

USING  DBUF ,4,5,6 
B  START  1 

DC  V(ICOP) 

CC  A ( ST  AP  T *4096 ) 

LA  B.INBUF  NUMPRI«8*(NPR-1>+JC0N 

MVC  TEMPI  3 )  ,.0<  B) 

MVC  TEMP2 ( 3 ( , 0 ( B  t 
L  9, TEMP 

SLL  9,  B 

SRL  9,16 

S  9, ONE 

SR  8,B 

M  8, EIGHT 

A  9 ,  JCON 

ST  9 , NUMPR  I 

L  9.NBET 

C  9 , NUMPP I 

BH  CDELTAR 

SPACE 

LA  B,toC273 

A  B.INDFX 

MVC  TEMP (3),0(6I 

L  9, TEMP 

N  9 ,  =X ’ FOCOOOOO ’ 

SRL  9,28 

ST  9,lC0nE  COMPUTE  THE  CODE  FCR  PRI 

C  9, THREE  SLAVEC  OR  NOT 

BE  CDELTAR 

C  9, SEVTN 

BE  CCELTAR 

C  9 ,  T  toO 

BE  CCELTAR 

L  9, TEMP 

N  9, =X 1 P  BCCOOOO ’ 

SRL  9,27 

ST  9.I27TB5  NBS  MCDE  INDICATOR 

L  9, TEMP 

N  9,  =  XT4CC0000» 

SRL  9,26 

ST  9,127386  ENDO-EXC  SCAN  INDICATOR 

L  9, TEMP 

N  9, =X ' 02000000 1 

SRL  9,25 

ST  9, 1 2  73B7  MBS  SCAN  MCDE  INDICATOR 
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GOOD  I 


S9 


SPACE 

L  9, TEMP 

N  9,=X»niCOCOOO' 

SRL  9 » 2  A 

ST  9.I273B8  DOUBLET  MODE  INDICATOR 

L  9, TEMP 

N  9  =X  *  PO IOOOOO  * 

SRL  9,20 

ST  9, I  278  12  NB/WB  INDICATOR 

SPACE 

LA  8,WC233  COMPUTE  GMT 

a  a.iNorx 

MVC  T EMP I  3 ) , 0 ( 8  I 

L  9, TEMP 

N  9  *  =X  *  I FCOOOOO  * 

SRL  9, 2A 

ST  9 , I H  R  S  STORE  HRS 

L  9, TEMP 

N  9, =X  ’PO3F00C0  * 

SR A  9,16 

ST  9 , I M IN  STORE  MINS 

L  9, TEMP 

N  9 , *X • POCO3F0G  * 

SR A  9,8 

ST  9,  I  SEC  STORE  SECS 

LA  8.WC23A 

A  8.IN0FX 

MVC  TEMP (11,0(81 

L  9, TEMP 

N  9, =X  •  7FEOCOGO  * 

SRL  9,21 

ST  9,  IMSFC  STORE  MSEC 

L  a, NSW 

C  B.ZERn 

RNE  S9 

SR  8,8 

SR  9,9 

LA  8.WC275 

A  8, INDEX 

MVC  TEMP ( 1 I , 0 ( 8 ) 

L  8, TEMP 

LA  9 , WC276 

A  9, INDEX 

MVC  TEMPI  3  1,0(9) 

L  9, TEMP 

SLL  9,1 

SLL  8,1 

SRL  a, 9 

SRDL  a, 6 

ST  9, FLMASK+A 

ST  8,FLMASK 

MV  I  FLMASK,X'A6' 

LD  0, FLMA  SK 

AD  0, D2EPG 

STD  0, XOPTIM 

SPACE 

LA  8 , WC26A  PRF  CALCULATION 
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» 


A 

8, INDEX 

NVC 

TEMP  ( 3 ) ,01 8 1 

L 

9, TEMP 

ST 

9 • WCRP6A 

LA 

8.WD273 

A 

8, INDEX 

NVC 

TEMPI  3 ) *0(6) 

L 

9, TEMP 

ST 

9,WCR073 

L 

9.WCRP6A 

N 

9 » =X • EFFF  EOOO  * 

SRL 

9,13 

KZSTMP 

ST 

9 , S  TEMP 

L 

9 ,  *F  •  10C0C000  * 

SR 

8,8 

C 

8 , STEPP 

ST 

SPACE 

9, STEPP 

L 

9,  INBl'F 

SRL 

9,31 

C 

9, ZERO 

BNE 

SPACE 

WB  AND 

L 

9 ,  WCRP73 

N 

9,=x‘nicooooo* 

SRL 

9 ,24 

C 

9, ZERO 

3E 

SPACE 

SL VCUP  1 

XOIV 

L 

B.FCUP 

XOIV1 

ST 

8,01  VSR 

B 

SPACE 

NEWPRF 

SLVDUB1 

L 

9,WCRP73 

N 

9 ,  -X  '  P8C0CC00 • 

SRL 

9,27 

C 

9, ZERO 

8E 

NBNWBN 

a 

XCIV 

K8NWBN 

L 

9 , WCR073 

N 

9 , -X • CO  1 OOCOO ' 

SRL 

9,20 

C 

9, ZERO 

BE 

NOD  I  VS 

L 

a, two 

B 

X0IV1 

KOCIVS 

L 

8, ONE 

B 

SPACE 

X0IV1 

WBANC 

L 

9.WCRP73 

N 

9  ,  *X ' 0 1 COOOOO • 

SRL 

9,2A 

C 

9, ZERO 

BNE 

SLVCUP2 

L 

8, TWO 

0 

X0IV1 

SLVDUB2 

L 

9, WCR073 

TRANSMITTED  PRF 

IK  NARROW  BANO 
BIT  8 

IK  DOUBLET  MOCE 

BIT  5 

IK  SLAVED  DOUBLET  MODE 
8 1  T  12 

NB/WB  E.O.P. 

NB  ONLY 

8 1  T  8 

IN  DOUBLET  NODE 
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NEWPRF 


NEXT* 


GOCON 


N 

9, *X  *08000000 

SRL 

9,27 

C 

9, ZERO 

BNE 

XDIV 

L 

8, TWO 

B 

XQIVi 

SPACE 

SR 

8,8 

L 

9.STEMP 

C 

8, DIVSR 

ST 

9, IPRF 

SPACE 

LA 

8,WC2-*7 

•4 

8, INDFX 

MVC 

TEMPI  3) ,0(0) 

L 

9. TEMP 

N  9 , *X  *  7FFFCOOO  * 

SRL  9,14 
ST  9.1AZ 

LA  8.WC276 

A  8 , INDFX 

"VC  TEMP ( 3 ) ,0 ( 8  ) 

L  9, TEMP 

N  9,*X«7FFFCCC0* 

SRL  9,14 
ST  9.IEL 

LA  8,WC26!» 

A  S.INOFX 

mvc  temp ( 7  j , o ( 8 ) 

L  9, TEMP 

N  9,  *X  *  FFFFEOOO  * 

SRL  9,13 

ST  9.TEMP2 

LA  8 , MC2F  7 

A  8.INOFX 

MVC  TEMPI 3 ) ,o ( 8  ) 

L  9, TEMP 

N  9,=X»FFFF0C00» 

SRL  9,16 

A  9, TEMP2 

S>  L  q,U 

S'  9, TEMP2 

LA  8 , WC2F6 

A  8, INDFX 

MVC  TEMP (3), 0(0) 

L  9, TEMP 

N  9, =x *  FFE  OOCOO  * 

SRL  9,21 

A  9, TEMP2 

ST  9, [RANGE 

LA  8 1 MC  1 1  5 

A  8, INDFX 

MVC  TEMP(3),0(aj 

L  9, TEMP 

N  9 , =X • POFFCOOO  * 

SR A  9,16 

ST  9, IPKPWR 


BIT  5 

IN  SLAVED  DGURLET  NCOE 
WB  ONLY 


STORE  A 2 


STORE  ELEV 


STORE  RANGE 


STORE  PEAK  POWER 
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LA 

8.MC269 

A 

8, INDEX 

MVC 

TEMP ( 3 ) , 0 1  8 ) 

L 

9 i TEMP 

C 

9,*F»n* 

BNL 

OOTG1 

N 

9,=X»7FFFFF00* 

SRA 

9,8 

LCR 

9,9 

B 

DOTG2 

COTGI 

SRA 

9,8 

CQTG2 

ST 

SPACE 

9 , I  ROOT 

LA 

8  ,  WC 1 1  7 

A 

8, INDEX 

MVC 

TEMPI  7 ) , 0( 8  ) 

L 

9, TEMP 

N 

9  ,  =  X  *  EFCOOCOO  * 

SRL 

9,24 

ST 

SPACE 

9, IMOVP 

L 

9, TEMP 

N 

9,=X*nCCCFF00» 

SRL 

9,8 

ST 

9,  IMOVQ 

SR 

9,9 

ST 

9,  ICFFST 

L 

9,1  COPE 

C 

9.THRFE 

BE 

CFFCO" 

C 

9, SEVFN 

BE 

OFFCOM 

B 

SPACE 

QFFSKP 

CFFCOM 

LA 

8  ,  WC278 

A 

8,  INOFX 

MVC 

TEMPI’), 018) 

SR 

9,9 

L 

9, TEMP 

C 

9, ZERO 

BNL 

RPLUS 

N 

9 ,  =  X  1  7FFFFF00 f 

SKA 

9,8 

LCR 

9,9 

B 

RNEG 

RPLUS 

SRA 

9,8 

RNEG 

ST 

SPACE 

9, ICFrST 

CFFSKP 

LA 

8 , WC26  3 

A 

8,  INDEX 

MVC 

TEMPI  1 ) , 0 1  8 ) 

L 

9, TEMP 

N 

9,=x'rocooooo* 

SRL 

9,26 

LA 

ll.PIFA 

LE 

0,0(9,11) 

STE 

0 , XPPAGC 

STORE  R-CCT 


ARE  PRIMARY  ANO  OFFSET  MOVING 


IS  OFFSET  hlNDOM  MOVING 


RANGE  OFFSET  FOR  SLAVED  WINDOW 


GET  VALLE  FROM  PIFA  TABLE 
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L 

9, TEMP 

N 

9  »  *  X • PFCCOOOO • 

SRL 

9,22 

LA 

1 1 , 0 1 F  A 

LE 

0,0(9,11) 

GET  VALUE  FROM  OIFA  TABLE 

STE 

O.XOPAGC 

L 

9,ZERr 

ST 

9,1 SWSSP 

ST 

9, ISWSSC 

ST 

9,  I  SSrRR 

LA 

8 , WD  239 

A 

8, INOFX 

MVC 

TEMPM)  ,0<  R) 

L 

9, TEMP 

f 

9,=X»OOCOC2O0' 

CHECK  BIT  23  (PFSA) 

C 

9, ZERO 

BE 

CKFSOP 

LE 

O.PFSA 

AE 

0, XPPAGC 

STE 

0,  XPPAGC 

ADD  IN  PFSA  VALUE 

CKFSCP 

L 

9,  TEM" 

N 

9, =X'COC00100* 

CHECK  BIT  24  (CFSA) 

C 

9, ZERO 

BE 

CKSSPP 

LE 

O.OFSA 

AE 

0,XCP«GC 

STE 

o.xcpagc 

ADD  IN  OFS A  VALUE 

CKSSPP 

L 

11, TEMP 

N 

ll,=X'OOdO2CO0' 

C 

1 1 «  *F  *  0  • 

ONE 

CK  S  SOP 

INOET' 

L 

8, ONE 

INDETERMINATE  SITUATION 

ST 

8,  ISSFRR 

B 

CCELTAR 

CKSSCP 

L 

11, TEMP 

N 

U,*X'0CAC100C' 

C 

1 1 ,  *F • 0  * 

BE 

INOET 

PPTEST 

LA 

9.HC239 

A 

9,  INOFX 

MVC 

TEMPI  7 ) , 0 ( 9 ) 

L 

10, TEMP 

AUX.MICR.WCRD  INTC  REG. 10 

LA 

9.WC2F2 

AUX.MICRCWAVE  WORD  INTC  REG 

A 

9, INOFX 

MVC 

TE  MP ( 3 ) , 0 ( 9 ) 

L 

11,  IE«P 

LA 

9,WC272 

A 

9,  INOFX 

MVC 

TEMP(3),0(9) 

RANGE  TR.WCRD  INTC  TEMP 

N 

10 , *X ' OC  802000 ' 

C 

1C , *X  *  00  800000 ' 

BNE 

S74 

LE 

O.PSSL 

ADO  IN  PSSL  (CCND.B) 

AE 

0, XPPAGC 

STE 

0, XPPAGC 

L 

9, ONE 

ST 

9,  ISWSSP 
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S7-. 

L 

8.NEWA 

UL3  OR  NEW  ATTEN. 

C 

8 » ZERC 

BE 

OPTEST 

L 

9, TEMP 

N 

9,=X‘rO080000* 

C 

9,=F'n* 

BE 

RDBKLC 

ATTENUATOR  READBACr 

N 

11,*X’OBOOOOOO» 

S7<  ARMED 

C 

lit  ZERO 

STATUS  READ  BACK 

BNE 

SLC 

NOATTLC 

LE 

O.PREVLC 

STE 

OtXPPAGC 

MVC 

JSWLCU)  ,ONE 

MVC 

I SSERR ( A  )  » ONE 

0 

OPTEST 

RDBKLC 

N* 

l lt^X'GAOOOOOC* 

S74  NOT  ARMED 

C 

1 1 1  ZEpO 

STATUS  READBACK 

BE 

NCATTLC 

B 

OPTEST 

SLC 

LE 

O.PSSA 

AE 

O.XPPAGC 

ADD  IN  PSSA  (COND.B) 

STCRLC 

STE 

OtXPPAGC 

MVC 

I  SWSSP ( 4  )  i CNE 

CPTEST 

LA 

9 ,  WC2  ^9 

A 

9, INDfX 

MVC 

TEMPI  1 ) tO( 9) 

L 

1C, TEMP 

AUX.MICR.WCRD  INTC  REG. 10 

LA 

A 

9,WC2‘i2 

9,  INDEX 

AUX. MICROWAVE  WORD  INTC  REG. 11 

MVC 

TEMP  I T ) , 0 1  9 ) 

L 

11, TEMP 

LA 

9.WC272 

A 

9,  INDEX 

MVC 

TEMPI 7  > ,019} 

RANGE  TR.WCRD  INTC  TEMP 

N 

10, =X'0CAC100C» 

C 

iC^X'OOACOOOO* 

BNE 

S75 

LE 

OtOSSL 

ACD  IN  CSSL  I COND.B ) 

AE 

OtXCPAGC 

STE 

O.XOPAGC 

L 

9, ONE 

ST 

9, ISWSSC 

S  75 

L 

8 , NEW A 

OLD  OR  NEW  ATTEN. 

C 

8 , ZERO 

BE 

OUT  1 

L 

9, TEMn 

N 

9,  =X ’OCCmCOCQ' 

C 

9 ,  *F • 0  ’ 

BE 

RD8KRC 

ATTENUATOR  READBACK 

N 

iitsX^OZOOOOOO' 

S 75  ARMED 

C 

11, ZEPO 

STATUS  READ  BACK 

BNE 

SRC 

NCATTRC 

LE 

OtPREVRC 

STE 

OtXCPAGC 

MVC 

JSWRC I  A ) , ONE 

MVC 

I SSERP I 4 ) ,CNE 

B 

OUT  1 
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WMtM  « 


RDBKRC 

N 

ll,=X*OIOOOOOO 

C 

lit  ZERO 

BE 

NCATTRC 

B 

CU1  l 

SRC 

LE 

O.OSSA 

AE 

0 • XCPAGC 

STCRCC 

STE 

0 « XCPAGC 

MVC 

I SWSSO ( 4 ) , ONE 

CUT  l 

L 

9, JSWLC 

C 

9, ZERO 

BNE 

OUT2 

LE 

O.XPPAGC 

SE 

0,=E* 16* 

STE 

O.XPPAGC 

STE 

0»  PREVLC 

CUT2 

L 

9, JSWRC 

C 

9, ZERO 

BNE 

ENOALPRT 

LE 

O.XOPAGC 

SE 

0,*E» 16* 

STE 

O.XOPAGC 

STE 

O.PREVRC 

ENCALER7 

MVC 

JSWLC ( 4  ) , ZERO 

MVC 

JSWHC ( 4 ) . ZERO 

L 

9  .  I  TB  AND 

C 

9f  ZERO 

BE 

NBANO 

LE 

2,RBIAS*16 

STE 

2,  TRBTAS 

L 

9.1  POLAR 

C 

9, ZERO 

BE 

lcpclar 

LE 

2.RBIAS+20 

AE 

2.TRBIAS 

STE 

2.TRBIAS 

L 

9,  ISWSSC 

C 

9, ONE 

BNE 

cceltar 

LE 

2. RBI AS*28 

AE 

2, TRBTAS 

STE 

2, TRBTAS 

B 

CCELTAR 

LCPOL AR 

L 

9, ISWSSP 

C 

9, ONE 

BNE 

cceltar 

LE 

2.RBIAS+24 

AE 

2, TRBTAS 

STE 

2, TRBTAS 

8 

cdeltar 

NBAND 

LE 

2, RBI  AS 

STE 

2, TRBTAS 

LA 

8.WD273 

A 

8, INOFX 

MVC 

TEMPI  31,018) 

L 

9, TEMP 

N 

9, =X ’ OCO ICOOO  * 

C 

9, ZERO 

S75  NOT  ARMED 
STATUS  REACBACK 

ADD  IN  CSSA  (COND.B) 


COMPUTE  RANGE  BIASES 
WIDE  BANC  TAPE 

OP  POLARIZATION 
ADD  WB  OP  BIAS 

ISWSSC  NAS  SET  IN  AGC  CCMP. 
*  1 » A  DC  32  DB  (OP) 

ADD  IN  CPS SA-  RBI  AS ( B ) 

ADD  IN  PSSA-RBIASI 7) 

NARROW  BAND 
CENTER  OR  EDGE  TRACK 
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BNE 

CKNBECGE 

B 

CKPCL AR 

CKNBECGE 

L 

a,  IRDCT 

C 

8, ZERO 

BH 

CKNBLOW 

LE 

2, RBI AS  +  A 

AE 

2,  TRBTAS 

STE 

2.TP8TIS 

fl 

CKPCL AR 

CKNBLCW 

LE 

2. RBI AS  +  8 

AE 

2, TRBTAS 

STE 

2, TRBTAS 

CKPULAR 

L 

?,  IPOLAR 

C 

9, ZERO 

BE 

CDEL  T  AR 

LE 

2*PB I  AS*  12 

AE 

2,  TRBTAS 

STE 

2,  TRBTAS 

CDEL  f  AR 

RETL 

TEMP 

CC 

F'  O' 

TEMP2 

CC 

F  •  0  • 

IXC 

CC 

F  •  0 ' 

NPTAPE 

CC 

F  •  0  • 

PKINUM 

CC 

F  •  0  • 

IPASS 

CC 

F '  0  • 

ISWSSO 

CC 

F  •  0  • 

ISWSSP 

CC 

F  •  0  * 

C  I  VSR 

CC 

F  •  0  • 

hORD6A 

CC 

F '  0 ' 

N0R07  3 

DC 

F '  0 ' 

STEMP 

CC 

F '  0 ' 

PREVLC 

CC 

E'O.O* 

PREVRC 

CC 

E'O.O' 

JShLC 

CC 

F '  0 ' 

JSWRC 

CC 

F  1  0 ' 

ZERO 

CC 

F'O' 

CNE 

CC 

F'  1* 

TWC 

CC 

F '  2 ' 

THREE 

CC 

F*  3' 

FOUR 

DC 

F '  A  • 

SEVEN 

CC 

F*  7' 

EIGHT 

CC 

F '  8 ' 

CIO 

CC 

F'  10' 

ClOO 

CC 

F'  100' 

CIOOO 

CC 

F'lCOO* 

CZERO 

CC 

co.o' 

FLMASK 

DC 

X' A60CCC000 

CBUF 

CSECT 

INBUF 

DS 

CL  3 

V.CI 

CS 

CL  3  PP 

CS 

CLA8 

V.CI8 

CS 

CL  3 

WD  19 

DS 

CL  3 

CS 

CL  2  7 

*D29 

DS 

<•3.3 

V.C30 

DS 

CL  3 

CS 

CL  8 1 

EDGE  TRACKING 
CENTER  TRACK 
CHECK  SIGN  CF  R  DC T 

LEADING  EDGE  BIAS 

TRAILING  EDGE  BIAS 
CHECK  POLARIZATION  DESIRED 

ACC  Ne  OP  BIAS 
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W058 

OS 

CL  17  1 

WD115 

OS 

CL  3 

WO  1 16 

CS 

CL  3 

WD1 17 

OS 

CL  3 

WD11B 

OS 

CL  171 

WO  1 75 

OS 

CL  1 7 1 

WD232 

CS 

CL  3 

WD233 

OS 

CL  3 

WD234 

OS 

CL  3 

CS 

CL  3 

W0236 

CS 

CL3 

WD237 

OS 

CL  3 

OS 

CL3 

W0239 

OS 

CL3 

WD240 

OS 

CL3 

WD241 

OS 

CL  3 

WD242 

OS 

CL  3 

OS 

CL  2  7 

W0252 

OS 

CL  3 

WD253 

OS 

CL  J 

CS 

CL  2  7 

WD263 

OS 

CL  3 

WD264 

OS 

CL  3 

WD265 

OS 

CL  3 

WD266 

OS 

CL3 

W0267 

OS 

CL  3 

WD  268 

OS 

CL  3 

WD269 

OS 

CL  3 

WD270 

OS 

CL  3 

WD271 

OS 

CL  3 

KD272 

OS 

CL  3 

WD273 

OS 

CL3 

WD27A 

OS 

CL  3 

WD275 

OS 

CL  3 

WD276 

OS 

CL  3 

WD277 

OS 

CL  3 

W027B 

OS 

CL3 

W0279 

OS 

CL  3 

WD280 

OS 

CL  3 

OS 

CL6369 

IA2 

OS 

IF 

I  EL 

OS 

IF 

INOEX 

DS 

IF 

IPPRCS 

OS 

IF 

IORS 

DS 

IF 

IRANGE 

DS 

IF 

IPKPWR 

OS 

IF 

IROQT 

OS 

IF 

I  ALT 

DS 

IF 

INDA2 

DS 

IF 

JNDAZ 

OS 

IF 

(NOEL 

DS 

IF 

IRB5A 

CS 

IF 

IRBB5 

OS 

IF 

IQPRCS 

OS 

IF 

I 2 A08  1 

OS 

IF 

I2A0B2 

CS 

IF 

I2A0B3 

DS 

IF 

I 24  IB  1 

OS 

IF 

I241B2 

OS 

IF 

PP  PHASE  0. 

CP  LOG  0. 

CP  PHASE  C. 
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I2AIB3 

OS 

IF 

XPPAGC 

DS 

IF 

!  BE  T  A 

DS 

IF 

NEW  A 

DS 

IF 

BAND 

DS 

IF 

NSW 

DS 

IF 

RBI  AS 

DS 

8F 

ISVPRI 

DS 

IF 

1  HR  S 

DS 

IF 

IMIN 

DS 

IF 

I  SEC 

DS 

IF 

IMSEC 

DS 

IF 

STAT 

DS 

21F 

TRBI AS 

DS 

IF 

ISTAT  1 

DS 

IF 

ISTAT2 

DS 

IF 

ISTAT3 

DS 

IF 

ISTAT* 

DS 

IF 

I  ALSW 

DS 

IF 

ISTSW 

DS 

IF 

NBWB 

DS 

IF 

ISIGNO 

DS 

IF 

127612 

DS 

IF 

JCCN 

DS 

IF 

NBEG 

DS 

IF 

NEND 

DS 

IF 

ITST 

DS 

IF 

NUMPRI 

DS 

IF 

XOPAGC 

DS 

IF 

IT6AND 

DS 

IF 

ITAPNC 

DS 

IF 

I  PRF 

DS 

IF 

I  POLAR 

DS 

F 

ISSERR 

DS 

F 

PIFA 

DS 

16F 

C  I  F  A 

DS 

16F 

PFSA 

DS 

IF 

CFSA 

DS 

IF 

PSS  A 

DS 

IF 

CSSA 

DS 

IF 

PSSL 

DS 

IF 

CSSL 

DS 

II F 

ICODE 

DS 

F 

127365 

DS 

F 

127366 

DS 

F 

I 27  36  7 

DS 

F 

I273B6 

DS 

F 

IMCVP 

DS 

F 

IMCVC 

DS 

F 

I OF F  ST 

DS 

F 

XOPT  I  M 

DS 

D 

IDAT 

DS 

682F 

END 
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APPENDIX  G 

SUBROUTINE  TIMDP  PROGRAM  LISTING 


SUBROUTINE  T I  MOP  I T I  ME , I  HR , MIN 1 1  SEC , I FR AC ) 
DOUBLE  PRECISION  T IME , T IME2 , XCON , FRAC 
DATA  XCGN/LOCOOOC.O/ 

I X=T IME 
IHR=I X/360C 
M I N=MOD ( IXf 360D)/60 
I  SEC-MOD ( I  X  »  60 ) 

T I ME2=0FLGAT ( IX)*XCCN 
T  ii  ME  =  T IME*XCCN 
FRAC=T IME-T IME? 

I FRAC=FRAC 

RETURN 

END 


APPENDIX  H 


SUBROUTINE  SMOOTH  PROGRAM  LISTING 


« 


r 


v 


SUBROUTINE  SMOOTH  ( N,  L  ,  X ,  NO ,  ZH  I 

C  SMOOTH  US/36C  CIMENSICN  RR.Cl.VEL  I NCREASEC  TO  700  3DJUN67 

C SMOOTH  OS/ 360  COMMON  CHANGED  ANO  ARGUMENTS  AS  HELL  12  JULY  66 

CSMCCTH  OS/ 360  20  JUNE  66 

CSMCCTH  CIMENStON  REOUCEC  TO  500  17  AUG  65 

CSMOCTH  26  MAY  65 


0  SMOOTH  REPLACES  THE  POSITION  MEASUREMENT  XIII  AT  THE  MIDPOINT  OF 

C  THE  INTERVAL  TU-NILESS  THAN  T  LESS  THAN  T(I*NI  BY  THE  CRCINATE 

C  OF  THE  PARABOLA  AT  THAT  POINT.  THE  SMCDTHEO  VALUE  XIII  THUS 

C  OBTAINEC  IS  EXPRESSED  AS  A  WEIGHTEO  AVERAGE  OF  THE  ME ASLREO 

C  VALUES.  THE  SLICING  PARABOLA  TECHNIOUE  CONSISTS  MERELY  OF 

C  APPLYING  THIS  PPERAIIGN  TO  EACH  INTERVAL  BETNEEN  TII-NIANO  T 1 1 *N» , 

C  SLICING  CNE  POINT  AT  A  TIME,  AS  I  *N* I ,N*2 , . . . ,L-N. 

C  ZH , THE  TIME  INTERVAL  BETWEEN  POINTS 

C  L.THE  TCTAL  NUMBER  OF  POINTS 

C  N.THE  NUMBER  OF  POINTS  CN  EACH  SICE  OF  XIII 

C  NO,  A  PARAMETER  WHICH  INDICATES  WHETHER  THAT  WHICH  IS  BEING  REAO 

C  IS  POSITION,  VELOCITY,  OR  ACCELERATION, 

f.  XII),  THE  X-COORDINATE 

C  YII),  THE  Y-COOROINATE 

C  Z ( I > ,  THE  Z-COORCINATE 

IMPLICIT  REAL*"  (A-H.O-XI 
COMMCN  RR 

DIMENSION  XI 1 , ,RRI I  I 
C 

M«2*N*1 

S*M 

JO«N*l 
KO«L-N 
SQ«S*S 
H  «  ZH 

5  IF  I  NO  I  5,10,20 

C  SET  UP  THE  COEFFICIENTS  FOR  POSITION 

5  0* 3.000/1  A. OCO*S* ( SC-5.0CO I 1 
01  -  SC-5. ODD 
C2  «  5 .OCD*S 
03  *  C 1*02 
05  >  3.DCD/C3 
02*D*0 

A  •  3.DDD*S0- 7.0C0 

8*0.0 

C--2D.D 

GC  TO  5757 

C  SET  UP  THE  COEFFICIENTS  FOR  VELOCITY 

10  0*12.0UD/IH*S*(SC-1.00C11 
01  *  SC-l.COC 
02  «  H*S 
03  *  0  I *C2 
05  *  I2.CCD/C3 
D2*C*0 
A*0.0 
B*I.O 
C*O.D 

GO  TO  5757 

C  SET  UP  T»  "OEFF I C I  ENTS  FOR  ACCELERATION 

20  D*3D.0C0/IH*H*S*(SC-l.>-.  MS0-5.0CC1) 

02*C*0 

A  •  l.DDO-SC 
9.0.0 
C* 12.0 

5757  DO  500  l*JC,KO 
SUMX«  A* X I  I  I 
00  600  X* 1 ,  N 
J*I*K 

LIK  *  I  -  X 
T*X 

V  *  -2.0DD»B*T 

SUMX*SUMX*(A*B*T*C*T*T 1*1X1 J)*X|LIXI )*Y*X(Ll»  I 
600  CONTINUE 

C  FORM  THE  SMOOTHEO  VALUE  FDR  EACH  COORDINATE 

ADJX*C*SUMX 
59  RR I  I  I  •  ACJX 
500  CONTINUE 
157  RETURN 
END 
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SUBROUTINE  DREFC  ( E t R* DEE « DRR )  VERSION1  6/15/70 

IMPLICIT  REAL  ♦  8  (A-H.O-Z) 

01  MENS  ION  CEllft,8),CR(l6,EI*ED(l6)»RDI8) 

DATA  DE/0.0  , U.O  *0.0  *0.0  ,0.C  *C.O  *0.0  .0.0  . 

10. 0  ,0.0  ,0.0  ,0.0  ,0.0  ,C.O  ,0.0  ,0.0  ,0.0313, 

20. 0303, 0.0292, 0.0287, 0.0282, 0.0272, 0.0262,0  0253,0.0263,0.0223, 

30. 0216, 0.0195,0. 0171, 0.0135, 0.0075, 0.0  ,0.0937,0.0868,0.0770, 

60. 0732, 0.0696,0. 0627, 0.0571, 0i0522, 0.068 0,0. 0612, 0.0385, 0.0337 , 
50.0278  0.0205,0.0105,0.0  ,0.1850,0.1520,0.1250,0.1160,0.1050, 

60.0906,0.0795,0.0708,0.0636,0.0523,0.0678,0.0605,0.0323,0.0229, 

70.0116,0.0  ,0.5310,0.3070,0.2120,0.1830,0. 1600,0.1280,0. 1060, 

80.0899, 0.0780,0.0612,0.0550 .0.0655,0.0356,0.0266,0.0120,0.0  , 

90. 7550,0. 3  720,0.2600,0.2020,0. 1 750, 0. 1 3 70  * 0 . 1  1 20,0.096 2*0. 08 1 1  , 

AO. 06 3 1,0. 0566, 0.0666,0. 036 1,0. 0250, 0.0 l 22, 0.0  ,0.9120,0.6110, 

BO. 2560, 0.2 160,0.1 860,0. 1620,0.1150,0.0967,0.0830,0.0663,0.0575, 

CO. 0672, 0.0 365, 0.0252, 0.0 122, 0.0  ,0.9700,0.6200,0.2600,0.2200, 

00.  1900, 0.  1660, 0.1  170, 0.0980, 0.08 60, 0.0653, 0.0586, 0.06 78,0. 03 69, 

EO. 0256, 0.0123,0.0  / 

DATA  CR/  C.O,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  C.O, 

1  0.0,  0.0,  P.C»  0.0»  0.0,  0.0,  C.O,  22.6,  21.5,  20.4,  19.9, 

2  19.6,  18.5,  17.6,  16.8,  16.1,  16.8,  16.2,  13.2,  12.0,  IC.6,  8.6, 

3  7.7,  67.3,  57.9,  50.2,  67.0,  66.1,  39.3,  35.6,  32.1,  29.3,  26.8, 

6  22.9,  19.7,  16.3,  12.7,  9.6,  8.1,132.0,  98.5,  77.6,  69.7,  63.2, 

5  52„9,  66.7,  3". 6,  33.6,  26.6,  23.9,  20.1,  16.6,  12.7,  9.6,  8.1, 

6360.0, 167.0, 103.0,  86.1,  73.6,  56.7,  66.2,  38.9,  33.6,  26.6,  26.0, 

7  20.2,  16.6,  12.8,  9.5,  8.2,605.0,170.0,106.0,  86.3,  73.6,  56.8, 

8  66.3,  38.9,  33.7,  26.5,  26.1,  20.3,  16.5.  12.8,  9.5,  8.2,621.0, 

9171.0,106.0,  86.6,  73.9,  57.1,  66.6,  35.0,  33.8,  26.8,  26.3,  20.5, 

A  16.6,  13. C,  n.8,  H. 6, 666.0, 172.0,105.0,  87.6,  76.0,  58.0,  66.6, 

8  39.2,  .36.0,  27. C,  26.6,  20.7,  16.7,  13.0,  10.0,  8.6/ 

DATA  EC, RTOto/O. 01, 2. 0,6. 0,5. 0,6.0, 8. 0,10. 0,12. 0,16.0, 18. ,20., 

126., 30. ,60. ,60. ,90. ,57. 29578/ 

DATA  RC/C.01, 10..3C. ,60.,200. , 600. , 1000 ., 2000. / 

IFIR.LE.0.01G0  TC  300 
RG*R/1. 8520*00 
00  100  I EO*  2 , 15 
I  -  17- 1  to 

IFIE.GE.EOI  I  1100  TC  12C 
100  CONTINUE 
1*1 

120  00  200  JRC=2,8 

J* 10—JK0 

IFIRG.Gt.ROI JIIGC  TC  220 
200  CONTINUE 
J*1 

220  IFIJ.EC.RIGO  TO  360 

ZR*CLOG(RG/RC( J) )/CLCG(RC(J*l l/RDI J 1 1 

1F(E.LE.C.0)G0  TC  320 

ZE*OLOG( E/EOI I ) ) /QLOG ( EO ( 1*1) /EDI  Ill 

0E1*((CE(  I*l,Jl-CE(I,jn6(l.-2R)6(0EII,J*n-CE(I,J))*2R16ZE 
0E2=( (0E( I , J*1 l-DEI l,J) )•(  l.-ZE)*(CEI I*l,J*l)-DE( I , J*1 1 l*ZE )*ZR 
'JEE  =  0E1*CE2*CE(  I ,  J  )  .  . 

PR1-M  (CR( I*1,J»-CR(I,J))*( l.-2Rl*(CR( I , J*1)-0R( I, J) >*ZR)*ZE 
OR 2* ( ( CR ( I , J* 1 1— CR ( I , J ) )*( I.-ZE )♦ COR  I  I *1 , J*l l-ORI I , J*l 1 16ZE)*ZR 
ORR= ( OR  »  *CR  2  *09 ( I , J)) 

GO  TC  600 
300  0EE*0. 0 

DRR-0.0 
GO  TO  60C 

320  DEE*OE(  I ,  J  )*U  Fi  ! ,  J*1)-DEI  I  ,  J  )  1*ZR 
ORR=OK ( I , J  )  ♦ ( C7 ( I , J*1 l-ORI I , J 1 )*ZR 

GO  TO  600 

3  60  OELT-IE-EDI  1)1/ (E0( 1*1 )-ED( I ) ) 

OEE*DEL ’•ICE  I  I ♦ 1, J l-OEl I , J ) 1 *DE ( I , J > 

ORR*CEL  T  * ( OR ( 1*1 ,J l-ORI I , J 1 1 ♦DR  1 1 , J 1 
600  ORR  ORR*. 3068D-03 
END 
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